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We propose a scheme for identifying an unknown Bell diagonal state. In our scheme the measurements are 
performed on the probe qubits instead of the Bell diagonal state. The distinguished advantage is that the quantum 
state of the evolved Bell diagonal state ensemble plus probe states will still collapse on the original Bell diagonal 
state ensemble after the measurement on probe states, i.e. our identification is quantum-state nondestructive. It 
is also shown finally how to realize our scheme in the framework of cavity electrodynamics. 
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I. INTRODUCTION 



X 



Entanglement is not only an essential feature of quantum 
mechanics, which distinguishes the quantum from classical 
world, but also a great resource in the fields of quantum infor- 
mation and quantum computation [1-3]. Particularly, entan- 
gled qubits prepared in the pure maximally entangled states, 
i.e., the Bell states, are required by many quantum informa- 
tion processes [4,5]. However, in the real world, a pure state 
in a quantum system will always evolves to a mixed one due to 
the unavoidable interactions with the environment. Thus, for 
practical purpose, applications of quantum information pro- 
cesses utilizing the mixed states are under consideration. 

Among the bipartite entangled mixed states, the Bell di- 
agonal state (BDS) plays an important role in quantum in- 
formation processing. It is widely used in the processes of 
quantum teleportation |^, quantum entanglement purification 
[7,8], quantum key distribution j^, etc. Moreover, the BDS 
is a simple but significant example in studying the nonclassi- 
cal correlation of a quantum mixed state [10-12], since there 
always exists a local transformation which can transform the 
given mixed state to a corresponding Bell diagonal form (v^. 
Therefore, identification of an unknown BDS is of great im- 
portance. Conventionally, the identification is achieved by 
the so-called state tomography technique [14-16], which per- 
forms the projection measurements on the unknown state di- 
rectly, and repeats the measurements on many copies of state. 
It is a drawback that after the projection measurements the 
state to be measured will collapse to one of the measurement 
basis. Thus the original state will be destroyed and become 
useless. 

Recently, schemes for quantum non-demolition measure- 
ment are proposed to detect unknown quantum state [17-19], 
by which the detected state will not be destroyed after the 
measurement. In this paper we present an alternative scheme 
for identifying an unknown BDS. In our scheme, we do not 
perform the projective measurements on the BDS directly, but 
on the probe qubits. According to the measurement outcomes 
of the probe qubits, we can acquire all the information of the 



unknown BDS. The distinguished advantage of our scheme is 
that the BDS is not destroyed by the measurements, since the 
evolved BDS plus the probe states will collapse to the origi- 
nal BDS after the measurement on the probe qubits. Contrast 
to the identification scheme with state tomography technique, 
which is achieved by sacrificing numerous copies of the un- 
known state, our scheme is economic and the resulting BDS 
is recyclable. The paper is organized as follows. In Sec. II, 
we explicitly demonstrate our scheme in theory. In Sec. Ill, 
we discuss the experimental realization of our scheme in the 
framework of cavity quantum electrodynamics (QED). The 
conclusion is drawn finally. 



'Electronic address: quaninfor mation@sina.com;ycs@dlut. edu.cn] 
^Electronic address: hssong@dlut.edu. cn| 



II. SCHEME FOR IDENTIFICATION OF UNKNOWN 
BELL DIAGONAL STATE 



In this section, we will illustrate our scheme explicitly. The 
BDS is a mixture of the well-known Bell states, it is param- 
eterized by four real numbers ci,C2,C3,C4 e [0,1], which 
satisfy the normalizing condition 2, c,- - 1. Generally, a 
BDS can be described as pn = 2/ c,|*P,)<*P,|, where = 

(|l>l|0)2 + |0>l|l>2)/V2,|>I'2) = (|l)l|0)2-|0)l|l>2)/^,|>I'3) = 

(|l>i|l>2 + |0>i|0>2)/ V2, and IT4) = (|l)i|l)2-|0)i|0)2)/ V2 are 
the Bell states, the subscript outside the ket denotes the label 
of qubit. Here |0) - [1,0]^ and |1> = [0, 1]^ are the computa- 
tional basis, the superscript T denotes the matrix transpose. 

In order to identify an unknown BDS, we need a probe 
qubit (labeled 3). The probe qubit interacts with the BDS and 
extracts the information from the BDS. We accomplish our 
scheme by three steps, in each step we can build an equation 
between the unknown parameters and the observable of the 
probe qubit. With the three steps done, we will obtain three 
independent linear equations, from which we can calculate the 
parameters. The three steps of our scheme are demonstrated 
in the following paragraphs. 

Step 1. Assume that the probe qubit is in state |0)3, thus 
the initial state of the joint system consists of the BDS and 
the probe qubit is given by - pn ® |0)3(0|. We perform an 
unitary operation f/', given as follows, on the joint three-qubit 
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As a result the state of the joint system evolves to p' = 
jyipOjyi 1 Qgjj obtain the reduced density matrix of the 
probe qubit by tracing over qubits 1 and 2 as foUows, 



1 _ I ci + C3 
Pi-\ 



C2 + C4 



(2) 



One can find that the information of the BDS is carried by the 
probe qubit. Performing a cr^ measurement on the probe qubit, 

we can obtain the following equation between the unknown 
parameters and the observable of the probe qubit, 



M' = Tr((T'^pl) - ci + C2 - C2 - C4. 



(3) 



Tracing over the probe qubit we can obtain the reduced den- 
sity matrix of the resulting BDS as follows. 



Pl2 = 5 



' ci+ C4 c\- c. 

C3 + C2 C3 - C2 

C3 - C2 C3 + C2 

c\- Cd, ci + C4 y 



(4) 



Note that underwent the f/^ operation, the resulting BDS be- 
comes different from the original state because the and 
IT3) ingredients have exchanged mutually. Fortunately, we 
can recover it to the original form by repeating the above- 
mentioned process once more with a new probe qubit to ex- 
change and 1^3) ingredients again. It is interesting that 
at the end of the recovering process, we can obtain the same 
reduced density matrix of the new probe qubit as shown in 
Eq. (2) and consequently yield equation Eq. (3) from the new 
resulting probe qubit. 

Step 2. In this step, the probe qubit is also initiahzed in |0)3, 
thus the joint system is in state pi2 ® |0)3{0|. We perform an 
unitary operation named on the joint three-qubit state, 
has the following form, 
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(5) 



After operation, the probe qubit evolves to the following 
form. 



2 _ ( Ci + C4 

C2 -H C3 



(6) 



Performing a cr^ measurement on the probe qubit, we can ob- 
tain the following equation, 



TrCoip,) = ci -H C4 - C2 - C3. 



(7) 



The resulting BDS underwent the f/^ operation is given as 
follows. 



' C3 + ci C3 - Ci 

C4 + C2 C4 - C2 

C4 - C2 C4 -H C2 

C3 - ci Ci+c\ ) 



(8) 



Similar to step 1 , we can transform the resulting BDS to the 
original form by performing on the joint system which is 
composed of the resulting BDS and a new probe qubit. Again 
the new resulting probe qubit ensemble will carry the infor- 
mation of the unknown parameters. 

Step 3. We perform an unitary operation on the joint 
BDS and probe qubit system, is given as follows. 
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Different from the previous two steps, here the probe qubit is 
initialized in the superposition state (|0)3 -1- 11)3)/ "V^. After 
operation, the probe qubit evolves to the following form. 



1 



1 



C\+C2- Ci - C4 



2 \ Ci -H C2 - C3 - C4 1 



(10) 



Now, performing a cr^ measurement on the probe qubit, we 
can obtain the following equation. 



= Tr(cr3P3) = c\ + C2 - - C4. 



(11) 



In the meantime, the resulting BDS has the following form. 



' C3 + C4 
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C4 - C3 



C\+C2 C\-C2 
C\- C2 C\ + C2 

V C4 - C3 C3+C4 ) 



(12) 



To recover this BDS back to the original state, we repeat 
on the joint system composed of the resulting BDS and a new 
probe qubit ensemble in state (|0)3 + 11)3)/ V2. Obviously, the 
new resulting probe qubit ensemble will also carry the infor- 
mation of the unknown BDS. 

Combining equations (3), (7), and (11), and taking into ac- 
count the normalizing condition 2, = 1> we can work out 
the parameters as: 



c\ = 



(13) 
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(16) 



Now we have succeeded in identifying an unknown BDS 
with the help of the probe qubit ensembles. Notably, due to 
the recovering process in each step, the final BDS is the same 
as the initial state. 

It is necessary give a discussion on the principles of our 
scheme. We emphasize that our scheme is based on the 
ensemble viewpoint, by which the BDS can be considered 
as a mixture of the four Bell states. The mixing propor- 
tion of each Bell state is denoted by the parameter c,. Each 
pair of qubits 1 and 2 fetching from the BDS ensemble will 
be randomly in one of the four Bell states. Without loss 
of generality, we take step 1 as an example to show how 
the probe qubit can extract information from the BDS en- 
semble. The expression of [/' can be rewritten as f/' = 

|lA2,()><iA2,()l + |lA2j)<<A2,ll + |lA4,o)<lA4,()l + |lA4J><<A4j|-#l,())<lA3J 

#i,i)<i/'3,ol - #3,o><iAi,il - #3,i><i/'i,ol, where = I^F/) ® 
One can find that if the state of qubits 1 and 2 is |*i'2) 
or |*1'4), it will remain unchanged and the probe qubit 3 will 
stay in |0)3; if the state of qubits 1 and 2 is I*!*:) (I^s)), it will 
change to (|*Pi)) and flip the probe qubit state from |0)3 
to 11)3. Repeatedly perform f/' on the joint three-qubit state 
by fetching new qubits from the BDS and the probe qubit 
ensemble, the resulting probe qubit ensemble will end in a 
mixed state ensemble which reveals the information of ci and 
C3 through the appearance probability of 11)3. To transform 
the resulting BDS back to the original form, we only need to 
repeat this process once more to make a simple exchange of 
I^Pi) and |*i'3). In steps 2 and 3, our scheme works similarly. 
As a consequence, the information of the unknown BDS is 
transferred to the probe ensembles. It is interesting that the 
resulting probe ensembles produced by the recover process 
are also useful. 

Let us look back to the expressions of U\ U^, and U^. 
These operators are essentially tripartite manipulations on 
qubits, and they can be formally factorized as U' = (C/J^ ® 
h){Ii ® U'2j){i =1,2, 3), where /i and I2 are the identity oper- 
ators of subsystems 1 and 2, respectively. The bipartite oper- 
ations f/'j3 and f/23 given as follows. 
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FIG. 1 : (Color online) (a) Schematic illustration for identification of 
unknown BDS based on the cavity QED system. The probe atom 3 
first interacts with atom 1 in cavity A and then is sent to cavity B 
to interacts with atom 2. Each atom in the cavity is driven by two 
classical lasers, (b) Atomic levels and transitions. 
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Based on these factorizations we can accomplish each step 
by sequentially performing bipartite manipulation on 
qubits 1 and 3, and f/^^ on qubits 2 and 3. That is to say we can 
perform only bipartite manipulations in the whole processing 
of our scheme, instead of tripartite manipulations which is dif- 
ficult to realize in experiments. The procedures are given as 
follows. Suppose that qubit 1 together with qubit 3 locates at 
place A, and qubit 3 locates at place B. In each step, we first 
perform operations C/j ^ on qubits 1 and 3, next send qubit 3 to 
place B, and then perform operations f/j^ on qubits 2 and 3. 
Finally, we perform measurements on the probe qubit. 



III. IDENTIFICATION OF BDS IN EXPERIMENTAL 
SCENARIO 

In this section we will discuss experimental realization of 
our scheme in the framework of cavity QED system. This 
experimental scenario is based on the case that two qubits 
are separated into different places, since two-qubit manipu- 
lation is more feasible than the three-qubit manipulation. The 
schematic illustration is shown in Fig. [Ua)- Assume that the 
BDS ensemble is shared by two participators, each of whom 
has an optical cavity A and B, respectively. The particles in 
the ensemble are considered to be three-level atoms with two 
ground states \a) and \b) and an excited state |e), see Fig. [Tfb)- 
The long-lived levels \a) and \b) represent states |0) and re- 
spectively. The probe atoms are identical to those in the BDS 
ensemble. The cavity couples the atomic transitions \a) <-> \e) 
and \b) «-> \e) with the coupling strength ga and gb, respec- 
tively. Additionally, two external driving lasers couple the 
transitions \a) <-> \e) and \b) <-> \e) with the Rabi frequencies 
and Q/,, respectively. All the atoms couple to the cavities 
via the same mechanism. 
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To start the scheme, we pick up a pair of entangled atoms 
from the BDS ensemble and put them into the corresponding 
cavities. The probe atom 3 is sent into cavity A firstly. The 
Hamiltonian in cavity A can be written as foUows, 

Ha = (Oe\e)j{e\ + (Oab\b)j(b\ + tOco'a 
J=13 

Hi^ae-'""' +gaa)\e)j{a\ 



+inbe-'"''' + gba)\e)j{b\, 



(20) 



where oj^ is the energy of |e) while cOai, is the energy of \b}, oja 
(ojh) is the frequency of the driver laser with Rabi frequency 
Cla i^b)< and a is the annihilation operator of the cavity. 

By setting 5\ - ojat - (pJa - <^h)l'^, we switch to an 
interaction picture with respect to Hq = Yji=\,3 <J^e\e)j{e\ + 
{coab - S\)\b)j{b\ + cOcO^a. Under the large detuning condition 
|5fl|,|5/,|,|Afl|,|A/,| » \ga\,\gb\,\^a\,\^b\, where 6a = w,. - w^, 

^b = COe-<^c-<^ab+Su = <^e-<^a, ^b = <^e-^b-<^ab+^U 

we can adiabatically eUminated the excited state \e)j [20-23]. 
If there are no photons in the cavity and the detunings satisfy 

\6a - ^b+ gllKl \h - K+ gll6a\ » \0.agblK\, \^bgal ^^b\ 

(we have assumed g\l6a = g^lSb), considering the subspace 
without real photons, we deduce the effective Hamiltonian as 

6a- ^b+ gll5a ^ 

where crt = \b)i{a\ and err = \a)i{b\, the coefficients are given 

as 



_ ga^b, 11. 



gh^a 1 1 
2 ^6b Aa 



(22) 



(23) 



B = 



^l^l/^^} . 1 



+ — ) + 

6a-6b'Aa Aj 6b- Aa+gl/Sa'Ab 6b 



1 



1 



( — + — ) 

6a- Aa+gll6a K 6a 
1 



(24) 



The effective magnetic field B can be tuned to be very close to 
zero by varying 6i. For ^b = gb^aij: + r)/ka(^ + ^)]. we 



^ * ^ "^[^ Oh" ^""^0 

can get the final effective Hamiltonian as follows. 



If we select the Rabi frequency as Clb = -gb^aij^ + 
^)/[^a(^ + ^)], we can obtain 



(26) 



where Ay = 2J\l{6a - Ab + gl/Sal At time t = ^^^,n = 
0, 1,2..., the time-evolution unitary operator e ''^n' coincides 
with the operator U^^. Then sent the probe atom into cavity 
B, we can realize the unitary operator U^^ by controlling the 
interaction time. 

In order to realize the operations U^^ and Ul^, we choose 
the laser frequencies as o)a = (Ob = oj. We switch to an inter- 
action picture with respect to Hq - 2j=i3 <jDe\e)j{e\ + (cDab - 
6\)\b)j{b\ + (jJcO^a, where the detuning 6\ is introduced to tune 
the effective magnetic field. Under the large detuning con- 
dition we can adiabatically eliminated the excited states. If 
there are no photons in the cavity and the detunings satisfy 
\6a -K+ gllSa\, \Sb - Afc + gll6a\ » \Q.agal K\, \^bgbl Ab\ (as- 
suming g^al^a - ghl^b), considering a subspace with no real 
photons we can obtain the following effective Hamiltonian, 



2(/i|a)i<a| + J2\b)x{b\){Ji\a)2{a\ + J2\b)2{b\) 



6a- K+ gl/Sa 



+B(cr, +o-\). 



(27) 



where cr^ - \b)i{b\ - \a)i{a\, and the coefficients are given as 
follows. 
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The effective magnetic field can be tuned to zero by varying 
6i. For Q.b = -^agaij^ + j;)/[gbij^ + ^)], we can obtain the 
following effective Hamiltonian, 



(31) 



(25) where = 2J\l{6a - A^ + gl/6a). It is obvious to see that the 
time-evolution unitary operator coincides with f/^j at 

time points t = ^^^i^^ (n = 0, 1, 2, ...). Thus we have realized 



where = 2J^/(6a - A;, -I- g^l6a). It is obvious to see that the 

unitary time-evolution operator e"''^"' is in accordance with the operation in cavity A, in the same way we can realize 
the unitary operator at t 



(2«+l);r 



0, 1,2..., thus we 



realize the unitary operation U\j,. To realize the operation f/js' 
we send the probe atom to cavity B, and drive atoms 2 and 3 
with the same lasers as done in cavity A. 



the operation in cavity B by sending the probe atom into 
cavity B. 

To confirm the validity of our approximation, we numeri- 
cally simulate the dynamics generated by the full Hamiltonian 
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FIG. 2: (Color online) Time evolutions of atom 3 calculated using 
full Hamiltonian (dashed blue line) and effective Hamiltonians (solid 
red line). The effective Hamiltonians for (a), (b), and (c) are 
and respectively. The initial states are chosen as 0.5(|a)i(a| + 
\b)m)® \ah{a\ for (a) and (b), and 0.25(|a)i<a| + \b)x{b\) ® (|a) + 
l^))3((fl| + (ft|) for (c). The parameters are chosen as = 1, f2„ = Sg,,, 
8a = 102g„, db = 122g„ A„ = 120g„, and A,, = 100g« for (a) and 
(b); g„ = 1, n„ = 5g„, <5„ = 102g„, Sh = 122g„, A„ = 100g„, and 
Aj = 120ga for (c). The inset shows the behaviors of the two curves 
in detail with time t 6 [500, 510] (in units of 1/ga). 

and compare it to the the dynamics generated by the effective 
Hamiltonian. In Fig. 2(a), we have plotted the time evolu- 
tion of (cr^). The numerical results show that the performance 
of (cr^) under the full Hamiltonian and that under H^^ agree 
with each other reasonably well. Similar agreement can also 
be seen in Fig. 2(b) and (c). Therefore, our effective model is 
valid. 

So far, we have realized all the unitary operations described 
by Eqs. (17)-(19) in the cavity QED system. Since the qubits 
are encoded in the ground atomic states and there is no real 
photons in the cavity, this experimental scenario is robust 
against the dissipative effects. 



Here we give a brief discussion on the experimental fea- 
sibility of the presented scheme. For an experimental im- 
plementation, the effective coupling strengths A^-, Ay, and A~ 
should be much larger than the cavity leaky rate Fc and the 
atomic spontaneous rate Fg. This requirements can be satis- 
fied in microcavities which have a small volume and thus a 
high quality factor Suitable candidates for the present pro- 
posal are, for example, the microtoroidal cavities which has 
cooperativity factor ^^/(FcF^) ~ 10^ and the ratio s/Te ~ 
10^ [23,24], where g is defined by g = max{g„,gh). Thus our 
scheme is feasible with cuiTent available systems. 



IV. CONCLUSION 

In conclusion, we have presented a scheme for nondestruc- 
tive identifying unknown BDS by measuring the probe qubits. 
This scheme is implemented in three steps. In each step we 
can build an equation between the unknown coefficients of 
the BDS and the observable of the probe qubit. Combining 
the three equations we can calculate the parameters. More- 
over, at the end of each step, the BDS ensemble remains in 
the initial state, therefore it is not polluted by the identification 
processing. We also consider the experimental realization of 
the scheme in the cavity QED system. By selecting appropri- 
ate Rabi frequencies of the driving lasers, we can realize the 
corresponding unitary operations, respectively. Our scheme is 
feasible with the current techniques. 
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